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ABSTRACT 

We present the results from coordinated X-ray observations of the ultraluminous X-ray source NGC 
5204 X-l performed by NuSTAR and XMM-Newton in early 2013. These observations provide the 
first detection of NGC 5204 X-l above 10 keV, extending the broadband coverage to 0.3 — 20 keV. 

The observations were carried out in two epochs separated by approximately 10 days, and showed 
little spectral variation, with an observed luminosity of Lx = (4.95 ± 0.11) x 10 39 ergs -1 . The 
broadband spectrum robustly confirms the presence of a clear spectral downturn above 10 keV seen 
in some previous observations. This cutoff is inconsistent with the standard low/hard state seen in 
Galactic black hole binaries, as would be expected from an intermediate mass black hole accreting 
at significantly sub-Eddington rates given the observed luminosity. The continuum is apparently 
dominated by two optically thick thermal-like components, potentially accompanied by a faint high 
energy tail. The broadband spectrum is likely associated with an accretion disk that differs from a 
standard Shakura & Sunyaev thin disk. 

Subject headings: Black hole physics - X-rays: binaries - X-rays: individual (NGC 5204 X-l) 


1. INTRODUCTION 

Ultraluminous X-ray Sources (ULXs) are off-nuclear 
point sources in nearby galaxies with observed X-ray lu¬ 
minosities Lx > 10 39 erg s -1 , exceeding the Edding¬ 
ton limit for a 10M 0 stellar-mass black hole (assuming 
isotropy). These high luminosities could be explained by 
a population of 10 2 — 10 5 Mq intermediate-mass black 
holes (IMB Hs) a ccreting at sub-Eddington rates (e.g. 
iMiller et, al.ll2004 1. Alternately, these luminosities could 
be explained by accretion onto stellar-remnant black ac- 
cretors (potentially reaching masses as large as ~1 00 M 0 ; 
iZampieri fe Robertsl 120091 Belczvnski et all [201011 if the 
emission is highly beamed fe.g. iKiii I EoQ 3), or if these 
sources ar e able to radiate in e xcess of their Eddington 
limit (e.g.lPoutan en et al.ll2007T). For r ecent reviews, see 
iRobertsI (120071 1 or lFeng fc Sorial (|201lH . 

ULX spectra in the soft X-ray band (< 10 keV) 
have been well-studied using telescopes such as XMM- 
Newton., Suzaku , and Chandra. Their spectral shapes 
appear to deviate substantially from those of known 
Galactic black hole binaries. A spectral turnover at 
<10 keV appears in most ULXs with sufficient sig nal-to- 
noise dStobbart et al.ll2o55 iGladstone et al.ll201lb . along 
with a soft excess usually modelled by a low tempera¬ 
ture < 0.3 keV blackbody disk component (IMiller et al.l 
I2004T) . This disk component, if produced by a stan¬ 
dard thin disk, could imply the presence of an IMBH 


accretor. However, the temperature-luminosity relation¬ 
ship of these sources does not appear to match the 
blackbody emissi on from standard accretion disks (e.g. 
iKaiava fo Poutanenl]2009i l. although the expected scal¬ 
ing may be partially recovered us ing a fixed absorp - 
tion column between observations (IMiller e t al.l l20 13aT) 
or us ing non-standard disk models dVierdavanti et ahl 
|200(J ). The low disk temperature can also be explained 
by a cool, optically th ick corona blocking t he inner 
disk from observation (jGladstone et al.l 12009 ). This 
corona would account for the continuum emission as 
well as the spectral turnover. Alternate possibilities are 
that the soft com pone nt orig inates from a strong out¬ 
flow (e.g. iPoutanen et al.ll2007? ) or blu rred line emission 
from highly ionized, fast-moving gas (|Goncalves fo Sorial 

EH). 

Until now, it has been difficult to distinguish between 
the spectral models due to the limited ~0.3-10.0keV 
bandpass over which ULXs have been studied. The dif- 
f erences become clearer with data above 10 keV (e.g. 

I Walton et al.l I2011H . a region of the spectrum that re¬ 
quires a focusing telescope with a broader bandpass. 
The Nuclear Spectro scopic Telescope Array ( NuSTAR ; 
lHarrison et al.ll2013fl . launched in June 2012, is the first 
orbiting telescope with hard X-ray focusing capabilities 
over a large 3 — 79 keV bandpass. With a similar effec¬ 
tive area to XMM-Newton at ~6keV, NuSTAR provides 
an ideal complement to the current soft X-ray observa- 
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TABLE 1 

Summary of X-ray data used in this analysis 


OBSID Detector Exposure (s) Counts 


Epoch 1 - 2013 April 19 


0693851401 

EPIC-pn 

13375 

9176 


EPIC-MOS1 

16396 

3010 


EPIC - MOS2 16458 

2846 

30002037002 

FPMA 

95964 

1871 


FPMB 

95799 

1907 


Epoch 2 — 

2013 April 29 


0693850701 

EPIC-pn 

10415 

6740 


EPIC-MOS1 

14036 

2129 


EPIC-MOS2 

14260 

2244 

30002037004 

FPMA 

88976 

1794 


FPMB 

88854 

1840 


tories for sensitive, broadband studies of ULXs. Indeed, 
over the past two years, NuSTAR , XMM-Newton, Chan¬ 
dra , Swift and Suzaku have undertaken joint observa¬ 
tions of several nearby ULXs (iBachetti et all 2013jj_ 2014 
iRana et aI1l2015t iWalton et al.ll2013bl . 12014 2015a lbf). 

This paper reports the results from observations of the 


ULX in NGC 52 04, a ne arby (d = 4.8 M pc) Magellanic- 
type galaxy ([Roberts &; Warwick! I2000D . NGC 5204 X-l 
has a typic al X -ray luminosity Lx ~ 2 - 6 x 10 39 erg s _1 


(Roberts et al. 

2004 

), and is well-studied below 10 keV 

(Roberts et al. 

2005 

. 2006) Vierdayanti et al. 2006j). It 


has been previously reported a s an IMBH candidate with 
long-term spectral variability dFeng fc Kaaretir2009D 

The paper is structured as follows: in section [2] we 
describe the observations and data reduction procedures, 
in section [3] we discuss the spectral analysis performed, 
and in section [I] we discuss the results and summarize 
our conclusions. 

2. OBSERVATIONS 

In April 2013 NuSTAR and XMM-Newton performed 
two coordinated observations of NGC 5204, approxi¬ 
mately 10 days apart. The NuSTAR exposures were 96 
ks and 89 ks, respectively, and the XMM-Newton expo¬ 
sures were 13 ks and 10 ks (EPIC-pn), and 16 and 14 ks 
(EPIC-MOS1/2). Details of the observations are sum¬ 
marized in Table [1] 

2.1. NuSTAR 

We reduced the NuSTAR data for each of the two focal 
plane modules (FPMA and FPMB) using standard meth¬ 
ods with version 1.1.1 of the NuSTAR Data Analysis 
Software (NuSTARDAS) and CALDB version 20130509. 
We ran the nupipeline tool to produce filtered event 
files, using all default options to remove passages through 
the South Atlantic Anomaly and periods of Earth occul- 
tation, and to clean the unfiltered event files with the 
standard depth correction, which substantially reduces 
the internal high energy background. We then extracted 
spectral products with nuproducts, using a 38" radius 
extraction region around the source, estimating the back¬ 
ground from a 113" radius region free of other point 
sources on the same detector as the target. The NuS¬ 
TAR data provide a reliable detection of NGC 5204 X-l 
up to ~20 keV. 


TABLE 2 

Parameters from a power law fit of both epochs 


Parameter Unit Epoch 1 Epoch 2 


XMM-Newton 


n H 

Wpi 

r 

X 2 /d.o.f. 

10 21 cm -2 
10“ 4 

0.47 ± 0.09 

3.36 ±0.12 
2.04 ± 0.04 
452/383 

0 49+ 0 10 
0.09 

3.34 ±0.15 
2.06 ± 0.04 
402/371 

XMM+NuSTAR 

n H 

N p i 

r 

X 2 /d.o.f. 

10 21 cm" 2 
10“ 4 

0.62 ± 0.09 

3.54 ±0.12 
2.12 ± 0.04 
591/451 

0 fit) - !" 010 
u ’ do -0.09 

o cc+0.16 

2.15 ±0.04 
538/436 


2.2. XMM-Newton 

The XMM-Newton data reduction was carried out 
with the XMM-Newton Science Analysis System (SAS 
vl2.0.1). To produce calibrated event files we used the 
tools epproc and emproc for the pn and MOS detec¬ 
tors, respectively. We then filtered out periods of high 
background according to the prescription in the SAS 
manualQ In evselect, we used the filters FLAG==0 && 
PATTERN<4 for EPIC-pn and FLAG==0 && PATTERN<12 
for the EPIC-MOS cameras. Spectra were extracted with 
evselect from a 30" radius region around X-l, and the 
background was estimated from a blank region of radius 
60" on the same detector, avoiding detector edges, bad 
pixels, and other visible sources. We also avoided the 
detector column passing through X-l, as recommended 
in the manual, to avoid the effects of charge spilling. An¬ 
cillary responses and redistribution matrices were gener- 
aged with arfgen and rmfgen, with the ELLBETA PSF 
correction enabled. 

3. SPECTRAL ANALYSIS 
3.1. General Procedure 

The spectral analysis for this work was conducted us¬ 
ing the Interactive Spectr al Interpretation System (ISIS) 
(iHouck fe Denicolal 120001) . ISIS was chosen over the 
more widely used XSPEC (lArnaudl [1996D for ease of 
programmability and its transparent use of parallelized 
fitting and error bar searches while also including all 
XSPEC models and tables. 


We modeled the neutral absorptio n co lumn usin g 


tbnew. a newer version of tbabs ( 

Wilms et al. 

20011, 

with the absorption cross-sections of 

Verner et al 

(199(|) 


and appropriate solar abundances. Cross-calibration be¬ 
tween the various detectors was addressed using a mul¬ 
tiplicative constant fixed to 1 for EPIC-pn and allowed 
to float otherwise; the calibrations of XMM-Newton and 
NuST A R are k nown to generally show a good agreement 
dMadsen et al.l I2015D . We performed fitting using \ 2 
minimization and quote errors as 90% confidence inter¬ 
vals unless stated otherwise. During our spectral analy¬ 
sis, all datasets were grouped to a minimum of 30 counts 
per bin to facilitate the use of \ 2 statistics. 

Fitting the XMM-Newton and broadband spectra in- 

1 http://xmm.esac.esa.int/ 
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Fig. 1.- Top panel-. Unfolded broadband spectrum of NGC 
5204 X-l. EPIC-pn is plotted in blue circles and EPIC-MOS1 and 
MOS2 are plotted in red stars and cyan squares, respectively. NuS- 
TAR FPMA and FPMB are shown in orange triangles and purple 
diamonds, respectively. Bottom panel'. Data/model ratio residu¬ 
als from an unabsorbed power law evaluated over the overlapping 
3.5 — 10 keV band, plotted over the full range. Data were rebinned 
to 100 cts/bin ( XMM-Newton) and 30 cts/bin ( NuSTAR ) for pre¬ 
sentation purposes. 

dependently with a simple power law indicates low vari¬ 
ability between the epochs, summarized in Tabic [2] The 
residuals for these fits behave very similarly in both 
epochs. For the remainder of the analysis we have there¬ 
fore combined the epochs using the HEASOFT tools 
addascaspec and addrmf to maximize source statistics. 

3.2. Modelling 

We first fit the combined spectrum restricted to the 
3.5 — 10 keV region of overlap between the XMM-Newton 
and NuSTAR using a simple power law continuum and 
Galactic ne utral absorp tio n co lumn, fixed at Nu = 1.66 x 
10 20 cm -2 (IKalberla et alJl2005f) . The result is a fit with 
xl = 1.18(199/168), and significant spectral curvature 
is evident when the model is evaluated over the broader 
0.3—20 keV spectral range, as shown in the lower panel of 
FigureH] The 0.3 — 3 keV spectral curvature may suggest 
an overall neutral absorption in excess of the Galactic 
column, so for the remainder of this work, we introduce 
a second absorption component intrinsic to NGC 5204. 
The column of this second neutral absorption model is 
allowed to vary. 

We fit the broadband 0.3 — 20 keV continuum using 
several models frequently used to describe ULX spec¬ 
tra. Initially, we examine six simple models: 1) a sim¬ 
ple power law; 2) a power law with exponential cutoff, 
XSPEC model cutof fpl; 3) a blackbody disk model with 
a radially variable temperatu re index, p diskpbb (a “slim 
disk” model with advection) (IMineshi ge et al.l ll994f): 4) a 
simple power law with an additional IShakura fe Sunvaevl 
(1973j)_ multicolor_blackbody disk component, diskbb 
( Mitsuda et al.lH984l) : 5) a power law with an exponen¬ 
tial cutoff and an additional multicolor blackbody disk; 
and 6) the same, replacing the exponential cutoff with a 


Energy (keV) 

Fig. 2.— Data/model ratios for some of the single-component 
models considered here: a simple power law (top), a power law 
with an exponential cutoff (middle), and an advection-dominated 
accretion disk with variable temperature profile (bottom). The 
color scheme is the same as Figure [l] 

broken power law, XSPEC model bknpower. Most mod¬ 
els we consider in this work are fit to both the broadband 
NuSTAR+XMM-Newton data and to the XMM-Newton 
data alone, for comparison, although we limit the major¬ 
ity of our description of the model fitting to the broad¬ 
band spectrum. 

Best-fit parameters for each model are summarized in 
Table [3j The simple power law model gives, as seen 
before, a poor fit with an “m”-shaped structure to the 
data/model residuals below ~10 keV and downward cur¬ 
vature in the 10 — 20 keV energy range. The power law 
with an exponential cutoff gives an improved fit with 
Ay 2 = 107 for one fewer degree of freedom (d.o.f.), 
but the “m”-shape is still visible in the residuals. We 
also attempted to model the spectrum using just an ab¬ 
sorbed multicolor blackbody disk, which gave a poor fit 
with y 2 = 6-18 (4582/741). At high accretion rates, 
t he expected emissi on likely deviates substantially from a 
IShakura fe Sunvaevl (|1973f) thi n disk, resulting in a sha l- 
lower temperature profile (e.g. lAbramowicz et ai1ll988l) . 
so we tried the diskpbb model, yielding a radial temper¬ 
ature profile p = 0.505 ± 0.004. While the diskpbb is 
a marked improvement over the simpler diskbb model, 
with y 2 = 1.21 (890/740), once again the “m”-shaped 
residuals remain, implying the need for two continuum 
components below lOkeV. The fits with these single com¬ 
ponent models (excluding the thin disk model, as the fit 
was very poor) are compared in Figure [2] 

The blackbody disk component introduced in the lat¬ 
ter three of our simple models considerably improves the 
low energy excess below ~ 2 keV. Adding the disk com¬ 
ponent to the powerlaw model provides an improvement 
of Ay 2 = 12 for 2 extra d.o.f., but the residuals for this 
model still show a downturn at high energies, so we also 
consider models including a disk component and a pow¬ 
erlaw with an exponential cutoff, and a broken power 
law. Both exponential cutoff and broken power laws pro- 
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vide a statistically good fit with y 2 = 1-04 (770/738 and 
771/737, respectively). The model including a blackbody 
disk and a broken powerlaw is perhaps a slightly bet¬ 
ter fit, with the data distributed more evenly about the 
model at high energies, although we cannot conclusively 
reject the cutoff powerlaw. The fits obtained with these 
two-component models are compared in Figured 

To explore more physical models, we investigate the 
possibility that the continuum is from cool, thin disk pho¬ 
tons Compton up-scattered in a hot corona. We model 
this component using the comptt model (jTitarchukl 
I1994T) . comptt is an analytic Comptonization model that 
assumes the seed photon spectrum follows a Wien law 
with some temperature Tq. Its use allows the temper¬ 
ature and the optical depth of the Comptonizing elec¬ 
trons to be fit as independent parameters. This model 
is frequently combined with a blackbody disk to repre¬ 
sent a standard disk-corona accretion geometry, with the 
Compton seed photon temperature linked to the inner 
disk temperature. Although we also present fits with 
the comptt model alone for completeness, we note that 
formally this describes a very extreme scenario in which 
essentially the entire X-ray emitting accretion disk is en¬ 
shrouded by the corona, which is likely unphysical. 

The fit with the comptt model provides a reasonable 
statistical fit with y 2 = 1-12 (824/739), but the “m”- 
shaped residuals seen previously are again apparent be¬ 
low lOkeV. Including a blackbody disk component gives 
a formally acceptable fit of xl = 1-07 (787/738), an im¬ 
provement of Ay 2 = 37 for one additional d.o.f., but 
evidence for an excess in the data over the model re¬ 
mains at high energies (see Figure 2] panel 2). This is 
also seen with si milar models in other NuSTAR UL X ob¬ 
servations (e.g. iWalton et al.ll2013bl 120141 1 201 511) . We 
address this excess with the addition o f a powe rlaw tail 
using the simpl convolution model dSteiner et al.ll2009i) . 
This powerlaw tail gives a fit improvement of Ay 2 = 
22 with two additional d.o.f. over the diskbb+comptt 
model, providing a statistically good fit with y 2 = 1.04 
(765/736) and resolving the high-energy excess seen pre¬ 
viously. In all of these models, the comptt component 
is cool and optically thick, resulting in a quasi-thermal 
blackbody-like spectrum. We note that the assumption 
of linking the Compton seed photon temperature to that 
of the observed disk emission may not be valid in the 
case of a central optically thick corona t hat obscures the 
inner disk (e.g. IGladstone et al.ll2009f) . However, we 
are unable to constrain these quantities independently if 
they are allowed to vary separately, and the fit remains 
unchanged, so we keep them linked for convenience. Re¬ 
placing comptt with a second diskbb component in this 
final model results in an equally good fit with y 2 = 1.04 
(766/737). Both models are illustrated in the bottom 
panels of Figure [I] A full list of best-fit parameters for 
these models is presented in Table 0] 

4. DISCUSSION AND CONCLUSIONS 

We have presented an analysis of the two coordinated 
NuSTAR and XMM-Newton observations of the ULX 
NGC 5204 X-l performed in 2013. The contribution of 
the NuSTAR data has allowed us to produce the first 
broadband spectrum of this source extending above 10 
keV. While NGC 5204 X-l is a source known previously 


TABLE 3 

Best fit spectral parameters for several common simple 

EMPIRICAL MODELS. 


Parameter 

Unit 

XMM-Newton 

XMM + NuSTAR 


tbnewx powerlaw 

UH 

Npl 

r 

X 2 /d.o.f. 

10 21 OH¬ 
IO - 4 

- 2 0.49 ±0.07 

3.50 ± 0.08 
2.05 ±0.03 
716/612 

0.65 ±0.06 
3.70 ±0.08 
2.14 ±0.02 
1030/741 


tbnew X cutoff pi 


n H 

10 21 cm 2 

0.49 ±0.07 

0 30+ 0 08 
U.OU—o 07 

N p i 

10~ 4 

3.51 ± 0.08 

o ^9+0.08 

°- u ^-0.07 

r 


2.05 ±0.03 

1.84 ±0.05 

E C ut 

keV 

500-392 

13 1+ 2 - 7 

J-O.i_i98 

X 2 /d.o.f. 


716/611 

923/740 


tbnew X di skpbb 


n H 

10 21 cm “ 2 

0.49 ±0.07 

0-37t°°e 

-^disk 

10“ 6 

4.08++ xl ° 5 

83.43t 2 ™ 2 

Tin 

keV 

7 99 + 6 . 8 I 

' 0.94 

3 C 5Q+ 0,31 
0.21 

V 


0.494 ± 0.004 

0.505 ± 0.004 

X 2 /d.o.f. 


705/611 

890/740 


tbnewX (diskbb+powerlaw) 


n H 

10 21 cm“ 2 

0.65t°;jt 

1 07+0.29 
1 - u '-0.26 

-^disk 


9.7S++ 9 

141 qn+ 0,32 
idi.yu_ 119 61 

Tin 

keV 

0.21 ±0.03 

0 1 oH-003 
U - 1Z —0.02 

N p i 

10~ 4 

9 04+°- 19 
z - y4 -0.20 

o 74 +0.l4 
'^-0.13 

r 


1 QQ+O-65 

t-oy_o.o6 

2.14 ±0.02 

X 2 /d.o.f. 


663/610 

1018/739 


tbnewX (diskbb+cutoffpl) 


n H 

10 21 cm 2 

0.28t°i« 

o.36t°;ji 

-^disk 


6-44 ±l;gg 

7-80± 4 ; 2 ° 

Tin 

keV 

0.28 ±0.03 

0.25 ±0.02 

Npi 

1 

O 

2-01 tg;£ 

9 qk?+0.26 

z.oU—027 

r 


0.68i°- 

1 OQ -1 " 016 
i - uy —0.18 

-Ecut 

keV 

q qi +1.78 
°- 0± — 0.96 

4 77 + 0 - 78 
' '-0.64 

X 2 / d.o.f. 


633/609 

770/738 


tbnew X (diskbb+bknpower) 



10 21 cm“ 2 

0.491°+ 

n 57+ 014 

U.O/-0.13 

-^disk 


6.52++ 

7 7z?+37.89 

1 • /o -3.40 

Tin 

keV 

0.25 ±0.03 

0.23 ±0.03 

N pl 

O 

1 

2 4i+ 0,29 
-0.33 

2.68+; 22 

Ti 


1 69 +0 10 

1,oy -0.13 

L80++ 

-^break 

keV 

5 07+ 0 - 72 
°- U '-0.81 

cr 01 +0.44 
°- ol -0.47 

r 2 


2 84 +0 ' 56 
z,o ^-0.45 

3 OO -1 " 016 

°- uu -0.15 

X 2 /d.o.f. 


633/608 

771/737 
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TABLE 4 

Best fit spectral parameters for several more physical 

MODELS. 


Parameter Unit XMM- Newton XMM + NuSTAR 


tbnew X comptt 


nu 

10 21 cm- 2 

o.oo + 008 

o.oo + 006 

Acomp 

1 

O 

3 5f)~*~ 105 

O.OU_o 35 

4 37+0-34 
-0.32 

To 

keV 

nil q+ 0.005 
U.llO-0 007 

0 H6+ 0 005 
u.iiq _ 0 006 

kT 

keV b 

O 9^+496.75 
o.z,o_o 73 

2 63” 1 ” 018 
z -°°-0.15 

Tp 


5-26 

6.01 ±0.27 

X 2 /d.o.f. 


680/610 

824/739 


tbnewX (diskbb+comptt) a 


n H 

10 21 cm" 2 

0.29l°;j® 

0-36±°;j| 

Adisk 


13 . 04 +IJg 

24.291^° 

T in 

keV 

0.23i°;« 

0.19 ±0.02 

Acomp 

1 

O 

2-89±° : « 

9 07+0.35 

z -°' - 0.30 

kT 

keV 


2 26”*~ 017 
z - zo - 0.14 

Tp 


9-33±m 

7 f) 7 + 0.50 
' ' — 0.46 

X 2 / d.o.f. 


633/609 

787/738 


tbnewX (diskbb+simpl*comptt) a 


nu 

10 21 cm -2 


0-30ln.i3 

Adisk 



9.86 lf 7 P 

Tin 

keV 


0.26 ±6.05 

Acomp 

1 

O 


3 14+0.90 

kT 

keV 


l-30t°j 2 

Tp 



10 . 32 ^ 2:19 

r c 



0 1 O+0.73 
0.10-2.08 

/scat d 



f) 61+0- 39 
u,D 1 -0.54 

X 2 /d.o.f. 



765/736 

tbnewx (diskbb+simpl*diskbb) 

nu 

10 21 cm -2 


0-25l°;j 2 

Adiski 



5.78t 2 ;° 5 ® 

T 

- L mi 

keV 


0.29 ±0.02 

Adisk2 

10“ 2 


2 24”*” b yb 
Z * Z4 -1.46 

T- 

- i m2 

keV 


1 13+0-45 
i-io_o .11 

r 



0 91 +0.31 
°' Z1 —0.73 

/scat d 



1.0— 0.70 

X 2 /d.o.f. 



766/737 


a The Comptonization input photon temperature has been set to 
the inner disk temperature. 

b Upper confidence bound at hard limit of kT = 500 for XMM- 
Newton data 

c Lower confidence bound at hard limit of T = 1.1 
d Upper confidence bound at hard limit of /scat = 1 



Energy (keV) 

Fig. 3. — Data/model ratios for the simple two-component con¬ 
tinuum models considered here: a simple powerlaw (top), a pow- 
erlaw with an exponential cutoff (middl e), and a broken pow erlaw 
(bottom), each combined with a lShakura fc Sunvaevl d 19731 ) thin 
disk model. The color scheme is the same as Figure [T1 
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Fig. 4. — Data/model ratios for the semi-physical models con¬ 
sidered: a single-component Comptonized continuum (top), thin 
a disk with an additional Comptonization continuum, the same 
model with an additional high-energy powerlaw tail, and the two- 
disk model with a high-energy power law tail (bottom). As before, 
the color scheme is the same as Figure [I] 


to display aperiodic spectral variability o n a time-scales 
of several days (e.g. iRoberts et al.l 120061) . we found low 
variability between the two observing epochs, separated 
by ~10 days, and therefore combined them to maximize 
count statistics for our spectral analysis. 

Prior work on NGC 5204 X-l using data from XMM- 
Newton suggested that it may be an IMBH of a few 
hundred solar masses, described by a cool blackbody 
dis k (kT ~ 0.2 keV) a nd a hard powerlaw tail (r ~ 
2; iFeng fc Kaaretl [200911 . Statistical evidence for cur¬ 
vature in the 2-10 keV bandpass has been seen in 


some previous observations (e.g. iStobbart et all 120061. 
which has generally been used to argue in favor of 
the high/super-Eddington interpretation, but not others 
(e.g. [Gladstone et al.ll2009i ). However, the limited band¬ 
pass of XMM-Newton meant even when this curvature 
was detected, it was not cle ar wheth er thi s represented 
a true spectral cutoff (e. g. iCaballero-Garcia fe; Fabiaiil 
120101 iWalton et al.l 120111) . Indeed, the limitations of 
a narrow bandpass are highlighted by the results pre¬ 
sented here, comparing the model fits to just the XMM- 
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TABLE 5 

Observed Fluxes for NGC 5204 X-l. 



Flux (erg s 1 cm 2 ) 

Luminosity (erg s 1 ) a 

0.3 - 10.0 keV 
10.0 - 20.0 keV 

0.3 - 20.0 keV 

(1.62 ±0.03) x lO" 12 
(l.T4±g-i|) x 10- 13 
(1.79 ±0.04) x lO' 12 

(4.47 ± 0.08) x 10 39 
(4.80±g;|f) x 10 38 
(4.95 ± 0.11) x 10 39 


a Observed luminosities (i.e. without any correction for neu¬ 
tral absorption app li ed) calculated for a distance of 4.8 Mpc 
IRoberts fc WarwicM2000f) . using the diskbb+simpl*diskbb model 


Newton data to those to the combined NuSTAR and 
XMM-Newton dataset. Far more models provide an ac¬ 
ceptable fit to the XMM-Newton data alone than to the 
broadband spectrum. 

The NuSTAR data presented here robustly demon¬ 
strate that the spectrum of NGC 5204 X-l displays signif¬ 
icant curvature above ~3keV, and is not powerlaw-like. 
This is similar to the re sults obtained in other NuSTAR 
ULX observations (e.g. Bachctti et al. 2013; Ran a et all 
120151 : IWalton et al.l l2013bL 12014 I2015allblh Flux calcu¬ 
lations further demonstrate that the proportion contri¬ 
bution of the hard X-ray emission from NGC 5204 X-l 
is relatively small, with only ~10% of the 0.3-20.0 keV 
flux emitted above lOkeV (see Table [5]). The ~5-6keV 
break argues against the presence of a ~ 10 3 M 0 IMBH 
accreting in the “low/hard state”, assuming such black 
holes would display similar accretion states to Galac¬ 
tic binaries, which would be ex pected appear powe rlaw- 
like i n the observed bandpass ([Remillard fe McClintock] 
I2006D . Using the best fit diskbb+simpl+comptt model, 
we calculate the observed 0.3 — 20 keV luminosity to be 
4.95 x 10 39 erg s _1 , (see Table [5]); the flux below lOkeV 
du ring this epoch is similar to that observed previously 
(|Roberts et al .1 [20041) . We conclude that NGC 5204 X-l 
is therefore likely a high-Eddington accretor with a more 
modest black hole mass. 

The broadband spectrum in this case is best fit by a 
three-component model, with two quasi-thermal compo¬ 
nents and a weak powerlaw-like excess at high energies 
(Figure [5t. We mod el the first component using a cool, 
IShakurafc Sunvaevi (119731 1 thin disk (Tj n ~ 0.25keV). 
The hotter components can be modelled as a cool, op¬ 
tically thick Comptonization model ( kT ~ 1.3 keV; r ~ 
10), dramatically different from the hot, optically thin 
coronae observed from sub-Eddington b l ack holes (e.g, 
Gieilinskijet^hlilOOOl Miller et alJl2013U: iTomsick et al.l 

2014t INatalucci et all 201 41 IBrennema.n et al.l I2014T) . or 
a second multicolor blackbody disk (Tj n ~ 1.15keV). 
Although poorly constrained, in both cases the best- 
fit photon index of the high-energy powerlaw tail is 
steep (r ~ 3), similar to the value of the photon 
index obsesrved in the steep pow er law state (T ~ 
2.5; iRemillard fc McClintock! 120061 ). In XSPEC syn¬ 
tax, the spectrum is best described by an absorbed 
diskbb+simpl*comptt or diskbb+simpl*diskbb model. 

In the context of high/super-Eddington accretion, a 
number of physical scenarios have been proposed for the 
emission components observed from ULXs below 10 keV. 
One such model invokes an optically thick Comptonizing 


corona that produces the 3-10 keV continuum and ob¬ 
scures the inner portion of the accretion disk, allowing a 
cool disk t emperature to be obse rved without requiring 
an IMBH ([Gladstone et al.ll2009l) . Our NuSTAR obser¬ 
vations likely do not favor this physical interpretation, 
as this model in turn seems to require a further Comp¬ 
tonizing region to explain the spectrum above 10 keV, 
which appears to have some similarity with the optically- 
thin Comptonization traditionally observed from Galac¬ 
tic binaries at high luminosities, calling into question 
the initial interpretation of the 3-10 keV cont inuum as 
the co rona. Furthermore, we note that iMiller et all 
(|2014T) argue that the parameters typically obtained with 
comptt fits to ULXs would imply very large size-scales 
for these coronae, and that it is difficult to envision a 
physical scenario that would result in such a large, uni¬ 
formly heated region . 

Alternately, iMiddleton et all (120111 ) have argued that 
the cool “disk” component could actually arise in a wind 
from a super-Eddington accretion disk, and that the 
hotter quasi-thermal model represents the spectrum of 
the inner disk itself. While we still lack unambigu- 
ous evidence of su c h wind s through absorption lines 
IWalton et al.l (120121 I2013all . this could be an effect of 
our viewing angle which may not intercept the winds 
if the they have a r oughly equato rial geometry (e.g. 
Middl eton et al.l [2015f) . In addition. iDexter fc Quataertl 
(1201211 have recently suggested that black hole binary 
accretion disks may exhibit significant inhomogeneities, 
resulting in the simultaneous presence of hot and cool 
regions within the same disk, and by extension an un¬ 
usual disk spectrum. Such “patchy” disk scenarios could 
arise as a natural signature of photon-bubble instabilities 
(|Gammiel 119981) proposed to trans port flux in a super- 
Eddington disk ([Miller et all 1 201 ill . Although there are 
differences in the detailed physics, both these models as¬ 
sociate the hotter thermal component with emission from 
the accretion disk, and the highest energy emission with 
optically thin Comptonization, and are consistent with 
our broadband observations. With current data, it is 
difficult to unambiguously associate model components 
with precise physical processes, particularly with regard 
to the soft thermal component; we cannot currently dis¬ 
tinguish between a disk or a wind origin here. Perform¬ 
ing broadband X-ray observations of NGC 5204 X-l at 
different flux levels to probe spectral variability and ex¬ 
amine how the different components evolve may be the 
key to distinguishing between these different models. 
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Fig. 5.— Left panel : the contribution of various spectral components for the diskbb+simpl*comptt model. The full model is given in 
black with the diskbb component shown in red and the comptt component (before modification by simpl) shown in blue. Right panel : 
same as the left panel but for the diskbb+simpl*diskbb model. The full model is again given in black with the first diskbb component 
shown in red and the second in blue. 
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